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Abstract 

Intermittency (externally induced) in the two-dimensional (2D) enstrophy cascade is shown to 
be able to maintain a finite enstrophy along with a vorticity conservation anomaly. Intermittency 
mechanisms of three-dimensional (3D) energy cascade and 2D enstrophy cascade in fully-developed 
turbulence (FDT) seem to have some universal features. The parabolic-profile approximation (PPA) 
for the singularity spectrum f(a) in multi- fractal model is used and extended to the appropriate 
microscale regimes to exhibit these features. The PPA is also shown to afford, unlike the generic 
multi-fractal model, an analytical calculation of probability distribution functions (PDF) of flow- 
variable gradients in these FDT cases and to describe intermittency corrections that complement 
those provided by the homogeneous fractal model. 
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1 Introduction 



Spatial intermittency is a common feature of fully developed turbulence (FDT) and implies that 
turbulence activity at small scales is not distributed uniformly throughout space. This leads to a 
violation of an assumption (Landau pQ ) in the Kolmogorov [2] theory that the statistical quantities 
show no dependence in the inertial range L S> £ S> r/ on the large scale L (where the external stirring 
mechanisms are influential) and the Kolmogorov microscale r\ = (z^/e) 1 / 4 (where the viscous effects 
become important), e being the mean energy dissipation rate. Spatial intermittency effects can be 
very conveniently imagined to be related to the fractal aspects of the geometry of FDT (Mandelbrot 
[3]). The energy dissipation field may then be assumed to be a multi- fractal (Parisi and Frisch [3J, 
Mandelbrot [5]). The latter idea has received experimental support (Meneveau and Sreenivasan 

In the multi-fractal model one stipulates that the fine-scale regime of FDT possesses a range 
of scaling exponents a E I = [a m in,«max]- Each a E I has the support set 8(a) C M 3 of fractal 
dimension (also called the singularity spectrum) f(a) such that, as i => 0, the velocity increment 
has the scaling behavior Sv(£) ~ £ a . The sets 8(a) are nested so that S(a') C 8(a) for a' < a. 

Experimental data on three-dimensional (3D) FDT (Meneveau and Sreenivasan [7j) suggested 
that the singularity spectrum function f(a) around its maximum may be expanded up to second 
order via the parabolic-profile approximation (PPA) 7 j|f - 

f(a)= f(a Q ) + l -f"{a ){a-a ) 2 (la) 

where, 

f(a ) = 3. (lb) 

On the other hand, spatial intermittency in two-dimensional (2D) FDT has been a controversial 
issue. Indeed, even the whole theory of 2D FDT (Kraichnan [10] and Batchelor [11]) had, until 
recently, remained almost an academic exercise, notwithstanding its possible connections with 
atmospheric and oceanic large-scale flows. 2D FDT has now been produced to a close approximation 

2 The PPA is equivalent to the log-normal model (Monin and Yaglom [S]). 

3 The PPA was shown to be a good fit for Lagrangian velocity statistics as well (Chevillard et al. [9]). 
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in a variety of laboratory experiments (Couder [12], Kellay et al. [13], Martin et al. [13], Rutgers 
[15] . Rivera et al. [16], Vorobieff et al. [17], Rivera et al. [18] and [E]). However, the Batchelor- 
Kraichnan theory for the enstrophy cascade in 2D FDT (with the energy spectrum E(k) ~ fc~ 3 jf| 
has essential differences with the Kolmogorov theory for the energy cascade in 3D FDT. The 
Batchelor-Kraichnan theory corresponds to flows with infinite enstrophy so there is no need for 
a nonlinearity-sustained cascade to transfer enstrophy across the inertial range to small scales to 
counter the dissipative action of viscosity there (Lopes Filho et. al. [H]). Consequently, the usual 
cascade mechanism is not operational in the Batchelor-Kraichnan theory for the enstrophy regime. 

Direct numerical simulations (DNS) of freely-decaying 2D FDT (McWilliams [22], Benzi et 
al. [23], Brachet et al. [24] and [25], Kida [26], Ohkitani [27], Schneider and Farge [28J) and 
forced-dissipative 2D FDT (Basdevant et al. [29], Legras et al. [30] and Tsang et al. [31]) showed 
intermittency caused by the presence of coherent structures. These structures inhibit the local 
inertial transfer of enstrophy via phase correlations across many length scales and produce energy 
spectra which are steeper than the Kraichnan-Batchelor spectrum E{k) ~ k~ 3 for the enstrophy 
cascade 1^| 

The theroretical issue of intermittency in the 2D enstrophy cascade in view of the regular 
behavior of 2D Navier-Stokes solutions is a delicate one. In addition to coherent structures, another 
cause of intermittency in the 2D enstrophy cascade can be contamination from the 3D effects in 
any real flow situation. In the geophysical context, it can be caused by the Ekman drag (|31j) 
which simulates the frictional planetary boundary layer (Pedlosky [33] ) . The Ekman drag acts also 
as a sink at low wave numbers to take out the 'condensate' from the inverse energy cascade. It 
is of interest to note that irrespective of the origin, an externally-induced intermittency is able 
to restore the usual nonlinearity-sustained cascade mechanism in the Batchelor-Kraichnan theory. 
Intermittency renders enstrophy of the flow finite so a nonlinearity - sustained cascade is very much 
needed to transfer enstrophy across the inertial range to small scales to counter the dissipative action 

4 The Batchelor-Kraichnan energy spectrum was recently shown (Eyink [20]) to correspond to the inviscid limit of 
the Leray solutions of 2D Navier-Stokes equations. These solutions, however, have vorticity fields which exist only as 
distributions. 

5 There is, therefore, some speculation whether different universality classes exist in the enstrophy cascade depend- 
ing on the particular initial conditions involved f|23|. |27] and |30|). On the other hand, the sensitivity of the 2D flow 
dynamics to the initial non-regularity of the flow raises questions about the whole universality concept in 2D FDT 
(Farge and Holschneider [32]). 
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of viscosity there. 

Paladin and Vulpiani [34] pointed out that the intermittency in the enstrophy cascade may 
be described in terms of a multi- fractal probability measure for the vorticity gradient^. The 
multi-fractal characterization of the intermittent enstrophy dissipation field was given by Mizutani 
and Nakano [57], Benzi and Scardovelli [38J, and Shivamoggi [39]. In the intermittent case, the 
vorticity at very small length scales (where viscous effects are strong and nonlinearities are weak) 
behaves like a passive scalar (Weiss |40j . Falkovich and Lebedev |41] and Nam et al. [42j ) and 
is advected by the large-scale flow structures |30| . This leads to thin sheets with large vorticity 
gradients (Saffman [35]) and selective rapid decay of vorticity in these layers because such regions 
experience typically stronger viscous diffusion than other regions. Conversely, vorticity-gradient 
layers (or divorticity sheets) are more likely to occur near vortex nulls. (This is very akin to 
the vortex-sheet formation near velocity nulls and current-sheet formation near magnetic nulls in 
magnetohydrodynamics (Shivamoggi et al. |44]-|46])). On the other hand, the consequent statistical 
dependence between the vorticity and the vorticity gradient at a point may be expected in turn 
to lead to non-gaussian statistics. Indeed, probability density functions (PDF) of enstrophy flux 
was measured in a freely-decaying 2D FDT by Kellay et al. [47] which was found to be highly 
non-gaussian. On the otherhand, DNS (Herring and McWilliams |48j . Borue [32], Maltrud and 
Vallis [50], Oetzel and Vallis [51]), soap- film experiments (Gharib and Derango [52], Kellay et al. 
|53j . Martin et al. [14] . Rutgers [E]), electromagnetically driven flow experiments (Paret al. |54J) 
showed restoration of the Kraichnan-Batchelor scaling behavior when the coherent structures are 
suppressed. 

The issue of whether intermittency (albeit externally-induced) in the 2D enstrophy cascade can 
maintain a finite enstrophy along with a vorticity conservation anomaly is of great interest - this is 
addressed in this paper. The possibility of common universal features (other than energy spectra) 
in the intermittency mechanisms of 2D and 3D FDT was raised by Dubrulle [55j and is also a topic 
of great interest. We explore here universal features in the intermittency mechanisms of 3D energy 
cascade and 2D enstrophy cascade (when it is intermittent via external induction). The PPA (1) is 



Indeed, experimental work of Jun and Wu [35] showed that large intermittency can be accounted for by the non- 
uniform distribution of saddle points in the flow (which are responsible for the energy transfer/dissipation (Daniel 
and Rutgers [36]). 
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used and extended to the appropriate microscale regimes for this purpose. We will also show that 
the PPA also affords, unlike the generic multi-fractal model, an analytical calculation of PDF's of 
flow-variable gradients in these FDT cases. 



2 3D FDT: Energy Cascade 

(i) Inertial Regime 

Let us briefly review and then extend the PPA applied to the inertial regime (Meneveau and 
Sreenivasan [7j, Benzi and Biferale |56j). 

According to the multi-fractal model for the pth order velocity structure function (Parisi and 
Frisch [4j), we have 

A P = (N P ) [ £ lpa+3 ~ na)] dfi(a) ~ (2a) 



where, 

4 1} = inf [pa + 3 - /(a)] (2b) 
and this minimum occurs for a = a* , which, according to the saddle-point method, is given by 

/'(«*) =P- (3) 

Writing (jlap and (jlb[) in the form - 

3 - /(a*) = a(a* - ao) 2 , a > (4) 

© yields, 

P 

a*(p)=a -—. (5) 
Using (@J) and ©, ([2a]) and (J2EJ) become 

4 1} =pao-^. (6) 
The parameter «q may now be determined (Benzi and Biferale [56]) by using the exact 3D 
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Navier-Stokes result (Kolmogorov [57]) - 

53 



4 1} = 1 (7) 



which, on application to ([6]), yields 



Using (|HJ), ([6]) becomes 



,p 3 ' 4a 
while © becomes 



13 

ao= 3 + 4a- (8) 



^ 1} = ? + ^(3-p)p (9) 



Mp) = 1 + ^M. do 

Next, using ([SJ and ©, (jl|) becomes 

2 

3-/(a,) = ^. (11a) 

Upon interpreting e as the energy transfer rate (in the inertial regime, as Kraichnan [58] pointed 
out, the energy transfer rate rather than the energy dissipation rate is the dynamically important 
parameter), we have from (9), 

(e 2 )~£~i;. (12) 
Comparing (12) with the log-normal result [8] - 

(e 2 > ~ (13) 



we have 



"=2a < 14 > 



Kolmogorov [59] and Obukhov [60J hypothesized that [i is a universal constant. 

It is of interest to note the sense in which the PPA (1) mimics a multi- fractal in describing 
the intermittency aspects of FDT. Comparing ([9]) with the multi-fractal result (Meneveau and 
Sreenivasan [6]) for the inertial regime - 

4 1) = 1 + ^(3-^(3-^/3) (15) 



we obtain 



Ava = 3(i-£). 



(16a) 



(|16a|) implies 



D = 3 



(16b) 



which is also confirmed by (11a) that yields 



/MO)) = 3 



(lib) 



/(a*(0)) being the fractal dimension of the support of the measure, namely, Do- It should be 
noted therefore that in the PPA the support of the measure is not a fractal; consequently, in the 
PPA the multi-fractality manifests itself via the way the measure is distributed rather than the 
geometrical properties like the support of the setQ In this sense the PPA is complementary to the 
homogeneous-fractal model (Frisch et al. [61J) in describing the intermittency aspects of FDT. This 
also implies that the failure to recognize the fractality of the support of the measure is apparently 
the cause of the well-known inability of PPA (See Castaing et al. |62j) to capture the quantitative 
aspects of intermittency adequately. 

(ii) Kolmogorov-Microscale Regime 

In order to determine universal features in the intermittency mechanisms of the various FDT 
cases, it is necessary to extend considerations to the Kolmogorov- microscale regime. Extension of 
the PPA to the latter regime has not been done. Let us now proceed to give this formulation. 

On extending the multi- fractal scaling to the Kolmogorov- rracrosca/e r/i, where, 



7 A multi-fractal generalizes, as Mandelbrot [5] clarified, the notion of self-similarity from sets to measures. 




(17) 



we have (Sreenivasan and Meneveau [63] and Nelkin |64j ) - 




(18) 
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where R\ is the Reynolds number - 



R 



Saddle-point evaluation of the integral in (18) yields 

(1 + a*)[p - /'(a*)] = Vol* - p + 3 - /(a*). (19) 

Using (HP, dE} leads to 

aal + 2aa* + (2p — 2aa>o — aafy = (20) 

from which, 



a,(p) = -l± 1 /(« + l) 2 --. (21) 

a 



Imposing the condition 



B ~ 1 (22) 

which, from (fl8]) . implies 

a*(0) = a (23) 
we see from (I2ip that the negative root needs to be discarded, and we obtain 



a*(p) = -l + A /(a + l) 2 --. (24) 

a 



On the other hand, using dU) and (p79|) . (fTHj) yields 



where, 



B p ~ fl^ (25) 



7 « = -[p + 2a{a*(p)-a }]. (26) 



In order to determine the parameter ao, the most pertinent framework for the Kolmogorov- 
microscale regime appears to be imposing the physical condition of inviscid dissipation of energy 
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(IDE) (Kolmogorov [57] )• This implies 



uB 2 ~ i?i 2 



constant 



(27) 



from which, 



75" 



1 = 0. 



(28) 



Using ([MD and QggJ, Q2SP yields 



1 3 

a °=3 + 4^ 



(29) 



which is identical to ([8|) that was obtained by imposing the exact 3D Navier-Stokes result ([7]) in 
the inertial regime! This is of course to be expected because the IDE is incorporated into the exact 
3D Navier-Stokes result (J7J). 
Using (|29D, (J23D yields 

I /A 3 \ 2 9n 

(30) 



while d2UJ then gives 



p + 2a\ 



16a + 9 
12a 



2p 16a + 9 
a 6 



(31) 



For large a, ([30]) and (H|) give the following asymptotic results 



(32) 



3-/K ) = g + o(^). (33) 

(|32p and ()33[) show that the zero-intermittency limit corresponds to a => oo, as before. 
(32) shows (as (10) does) that 

a* < ~, Vp > 2, (34) 

implying of course the strengthening of the velocity-field singularities by intermittency in the 
Kolmogorov- microscale regime! 
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(iii) Probability Distribution Function for the Velocity Gradient 

The multi-fractal model is known not to afford an analytical calculation of PDF of velocity 
gradient in intermittent FDT (Benzi et al. [65]). We now wish to show that the PPA is fruitful 
on this aspect. The physical principle underlying the calculation of the intermittency correction 
to the PDF of velocity gradient in the PPA turns out to be however the same as the one (namely, 
IDE) underlying the homogeneous-fractal model used in [65J . 

Noting the scaling behavior of the velocity gradient (Frisch and She |66j ) - 



(35) 



vo being the velocity increment characterizing large scales, and assuming vq to be gaussian dis- 
tributed, i.e., ,,2 

i> )~e ^} (36) 



we observe 



(37) 



So, a*{p) corresponds to a*(p) where p is the solution of 



p = l + a*(p). 



(38) 



Using (|38|) . and assuming a to be large to simplify the calculations, we have from 



M P ) = - 3 - Ta +o 



(39) 



Using (|39|) . the PDF of the velocity gradient [66] 



[i-"«(p)] 



^{l-q«(p)}| 3 |{l+q«(p)} 



(40) 



becomes 



P(s,a*(p))~ ( ^ 



>3 T 8o, 



„(§+*«) i s |U - *s) 



(41) 
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Incidentally, using (|39|) . (|38|) gives 



P=t-^ + o(\) (42) 



3 4a \a 2 

which is of course the exponent of |s| in the argument of the exponential in (|4ip as to be expected 
from (37). Note the accentuation of the non-gaussianity of the PDF due to intermittency, as 
also indicated by the homogeneous-fractal model [65] which, as pointed out before, is however 
complementary to the PPA in describing the intermittency aspects of FDT. 

3 2D FDT: Enstrophy Cascade 

The evolution of vorticity in a 2D fluid flow is governed by 

^ + (v • V)w = vV 2 uj (43a) 
ot 

and is based on the competition between the viscous diffusion and the advection processes. 
Writing (43a) in the form 

^ + V • (vu>) = vV 2 u (43b) 

we observe that vorticity is globally conserved in the absence of viscous diffusion, so the local 
maxima of vorticity can grow and possibilities of generation of non-Gaussian statistics for the 
vorticity PDF exist. 

On the other hand, taking the gradient of (43a), we obtain 



^(Vw) + (v • V)(Va;) = - (Vtfe|£ + + ^V 2 (V^). (44) 

(44) shows that the vorticity gradient is not conserved along the Lagrangian trajectories even 
in the absense of viscous diffusion, so the local maxima of the vorticity gradient field can grow. 
This would imply that the statistics of the vorticity gradient can also become non-Gaussian. 

We therefore consider intermittency (albeit externally-induced) in the enstrophy cascade of 2D 
FDT and extend the PPA to formulate this. 
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(i) Inertial Regime 

The multi- fractal model for the pth. order velocity structure function gives (Shivamoggi [39] ) 

A p ~ £& } (45) 

where, 

4 2) = inf \pa + 2 - /(a)] (46) 
and this minimum a = a*, according to the saddle-point method, corresponds to 

/'(«*) =P- (47) 
Assuming a PPA for the 2D case, we have - 

2 -/(a,) = a (a*-a ) 2 , a>0 (48) 

Using (48), (47) yields, 

a*(p) = "0 - 7T ■ (49) 
za 

Using (gSD and flU, ([Ml) becomes 

4 2) = - C ■ (50) 

Using the exact result for the enstrophy cascade in 2D FDT (Eyink [67] and Lindborg [68j]f - 

= 3 (si) 

the parameter ao may be determined - 

a = 1 + -f- ■ (52) 
4a 



Equation (51) is also predicted by the multi-fractal model for the enstrophy cascade in 2D FDT (Shivamoggi [39] 
and [69]) and is in agreement with the rigorous inequalities established by Eyink '70 and has been verified by DNS 
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Using (55), (53) becomes 



while (52) becomes 



4 2) =P+^(3-p)p (53) 



aM = l + ^j^- (54) 

(54) shows that 

a«,(p)<l, Vp>2 (55) 

implying of course the strengthening of the vorticity-field singularities by intermittency! Indeed, 
(53) predicts an energy spectrum - 

E{k) ~ £T 3 ~^ (56) 

which is steeper than k~ 3 , as required. 

Observe that the intermittency corrections in (50) and (51) are identical to those for the energy 
cascade in 3D FDT, namely, (8) and (9)! 

On the other hand, (53) implies for the structure function of the vorticity field, the scaling 
behavior - 

(|<^| p ) ~ 6^ (57) 



where, 



<« = -g(p-3). (58) 



Comparing (57) with the multi- fractal result (Benzi and Scardovelli [38j)| 

(P = -\( P -3)(2-D p/3 ) (59) 



9 (58) is equivalent to [39] 

& } = p-i(p-3)(2-Z3 p/3 ) 

which implies that the energy spectrum in the enstrophy cascade cannot be steeper than fc -11 / 3 (corresponding to 
Dp/3 = 0) in agreement with Sulem and Frisch [71] and Pouquet [72]. Further, the result E(k) ~ fc _11//3 agrees with 
the result of Gilbert [73; which considered the dissipative structures to be line elements (with zero fractal dimension 
for their projections on the plane) - these lines are centers of accumulation of singularities associated with spiral 
vortex sheets in Gilbert's model. 
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we obtain 




(60a) 



Observe that the intermittency correction in (60a) is again identical to that for the energy cascade in 
3D FDT, namely, (16a)! Incidentally, the necessity of negative Holder singularities of the vorticity 
field to produce intermittency in the enstrophy cascade and to preclude enstrophy conservation in 
the inviscid limit was noted by Eyink |67j . 
On the other hand, (60a) implies 



which signifies that in the PPA the support of the enstrophy-dissipation field in 2D FDT is not a 
fractal. This has been confirmed for 2D FDT however by the DNS [30]. 

These results appear to signify universal features in the intermittency mechanisms of 3D energy 
cascade and 2D enstrophy cascade, 
(ii) Kraichnan Microscale Regime 

On extending the multi- fractal scaling to the Kraichnan microscale rj2 (Shivamoggi |74j). where, 



D 



= 2 



(60b) 




(61) 



t being the mean enstrophy dissipation rate, we have [39] 




(62) 



where Ri is the Reynolds number for the 2D FDT - 



R2 



( r L 6 )V3 



Saddle-point evaluation of the integral in ([621) yields 



(1 + a*)[p - /'(a*)] = pa* - 2p + 2 - /(a*). 



(63) 
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Using (gHD, (JB3J leads to 



aal + 2aa* + (3p — 2aao — aafy = (64) 



from which, 



Imposing the condition 



on (62), we have, 



a,(p) = -l± A /(a + l) 2 - — • (65) 

a 



C ~ 1 (66) 



a*(0) = ao- (67) 
Using (67), we see that the negative root in (65) is to be discarded, and we obtain 



n,(/;) = -l+ • (no • 1)-' ^. (68) 



On the other hand, using (48) and (63), (67) yields 



(2) 

C p ~ i?^ p (69) 



where, 

7 ( 2 ) = -[p + 2a{a*(p)-a }]. (70) 

In order to determine the parameter ao, the most pertinent framework for the Kraichnan- 
microscale regime in the 2D enstrophy cascade appears to be imposing the physical condition of 
inviscid dissipation of enstrophy (IDE) [39J This implies 

(2)i 

vC 2 ~ Rl 2 ~ const (71) 



10 Indeed, Polyakov [75] argued that the existence of the enstrophy cascade is predicated on the enstrophy conser- 
vation law anomaly, while recent numerical computations (Dmitruk and Montgomery [76]) showed that the tendency, 
if any, of enstrophy dissipation to go to zero in the inviscid limit is a very weak one. 
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from which, 



Using (68) and (70), (72) yields 



1? ] "1 = 0. (72) 



«o = l + ^ (73) 



which is identical to (52) that was obtained by imposing the 2D multi-fractal model result (51) in 
the inertial regime! This appears to indicate that the IDE has been incorporated into the result 
(51) in a manner similar to the case with the Kolmogorov exact result in the 3D incompressible 
case. Thus, intermittency (albeit externally-induced) can maintain a finite enstrophy along with a 
vorticity conservation anomaly without contradicting the rigorous result in [20] and [21]. 
Using (73), (68) yields 

<fcW = -l + ^( 2+ ») , -f (74) 

while (70) then gives 

'3 + 8a\ 2 3p 3 + 8a 



^( 2 ) = - 



V + 2a<i 



4a J a 



(75) 



The intermittency corrections Aa*(p) to the scaling exponent for the 3D (given by (30)) 

and the 2D (given by (74)) are sketched in Figures la and lb for different values of the intermittency 
parameter a. Observe that the velocity-field singularities are strengthened in the microscale regimes 
by the intermittency effects. Intermittency corrections Aa*(p) for 2D cases are smaller than those 
for 3D cases. Further, in the weak- intermittency limit (a large), the intermittency corrections 
Aa*(p) for the 3D and 2D cases are almost identical. 

For large a, (74) and (48) give the following asymptotic results - 

«.(,,)= l + (l- p )l +0 (i) (76) 



Observe that the intermittency corrections in (76) and (77) for the 2D enstrophy cascade microscale 
regime are identical to those in (32) and (33) for the 3D energy cascade microscale regime! This 
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appears to confirm further the universal features in the intermittency mechanisms of 3D energy 
cascade and 2D enstrophy cascade pointed out earlier. 



(iii) Probability Density Function for the Vorticity-Gradient 

Noting the scaling behavior of the vorticity gradient (Shivamoggi [39] 



r ~ — 7j ~ v n +a * v 1+a 



(78) 



and assuming ([36]) again, we observe 



(79) 



So, a*(p) corresponds to where p is now the solution of 



P = gl 1 + «* 



(80) 



Using (80), and assuming again a to be large, we have from (76), 



«.(p) = 1 - - + o , „ 



4a V a' 



(81) 



Using (81), the PDF of the vorticity gradient [39] 



12-"* (a: 

P(r,a*(p)) ~ ( j-T I ( 



[/ §{2-c,(p)} | ,§{l+c,(p)} 



(82) 



becomes 



// \ 3 ( 1+ 4a ) 

P(r, ~ ( j-r I ' 



I/ §( 1 +^-)| r ,|§( 2 -ir) 



Using (81), (80) gives 



(83) 



(84) 



3 V 4a 

which is again the exponent of |r| in the argument of the exponential in (83), as to be expected 
from (79). Note again the accentuation of the non-gaussianity of the PDF due to intermittency, as 
also indicated by the homogeneous-fractal model [74J which is however complementary to the PPA 
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in describing the intermittency aspects of FDT. 

4 Conclusions 

The theoretical issue of intermittency in the 2D enstrophy cascade in view of the regular behavior 
of 2D Navier-Stokes solutions is a delicate one. It is of interest to note that, irrespective of the ori- 
gin, an externally-induced intermittency is able to restore the usual nonlinearity-sustained cascade 
mechanism in the Batchelor-Kraichnan theory. Intermittency renders enstrophy of the flow finite 
so a nonlinearity sustained cascade is very much needed to transfer enstrophy across the inertial 
range to small scales to counter the dissipative action of viscosity there. Further, as we have shown 
in Section 3 (ii), intermittency can maintain a finite enstrophy along with a vorticity conservation 
anomaly without contradicting the rigorous result in [20] and [21]. 

Intermittency mechanisms of 3D energy cascade and 2D enstrophy cascade appear to have 
certain universal features notwithstanding very different physics underlying these FDT cases. This 
probably results, following Kadanoff 's [77] speculation, because of consideration of the limit R => oo, 
which happens to be the critical point for FDT. Indeed, universal features for very diverse FDT cases 
become apparent at the critical point (Shivamoggi [75]). On the other hand, via (14) this appears 
to support the universality of the log-normal exponent \i hypothesized by Kolmogorov [59J and 
Obukhov [60]. It remains to be pointed out, however, though the PPA appears to have the capacity 
to provide considerable insight into the qualitative aspects of intermittent FDT quantitative aspects 
have been a different story for the 3D FDT problem (see Castaing et al. [62]). This is probably 
traceable to the failure of PPA to recognize the fractality of the support of the measure, as shown 
in this paper. 

Acknowlegements 

This work was carried out during the author's visit to the International Centre for Theoretical 
Physics, Trieste, Italy. The author is thankful to Professor Katepalli Sreenivasan for his valuable 
remarks and discussions. 



18 



References 

[1] L. D. Landau and E. M. Lifshitz: Fluid Mechanics, II Ed., Pergamon Press (1987). 
[2] A. N. Kolmogorov: Dokl. Akad. Nauk. SSSR 31, 19 (1941). 

[3] B. Mandelbrot: in Turbulence and Navier-Stokes Equations, Ed. R. Temam, Lecture Notes in 
Mathematics, Vol. 565, Springer- Verlag (1975). 

[4] G. Parisi and U. Frisch: in Turbulence and Predictability in Geophysical Fluid Dynamics and 
Climatic Dynamics, Ed. M. Ghil, R. Benzi and G. Parisi, North-Holland (1985). 

[5] B. Mandelbrot: in Fractals in Geophysics, Eds. C. H. Scholz and B. Mandelbrot, Birkhauser 
(1989). 

[6] C. Meneveau and K. R. Sreenivasan: J. Fluid Mech. 224, 429 (1991). 

[7] C. Meneveau and K. R. Sreenivasan: Nucl. Phys. B Proc. Suppl. 2, 49 (1989). 

[8] A. S. Monin and A. M. Yaglom: Statistical Fluid Mechanics, Vol. 2, MIT Press (1975). 

[9] L. Chevillard, S.G. Roux, E. Leveque, N. Mordant, J.F. Pinton and A. Arneodo, Phys. Rev. 
Lett. 91, 214502 (2003). 

[10] R. H. Kraichnan: Phys. Fluids 10, 1417 (1967). 

[11] G. K. Batchelor: Phys. Fluids 12, Suppl. 11-233 (1969). 

[12] Y. Couder: J. Phys. Lett. 45, L-353 (1984). 

[13] H. Kellay, X. L. Wu and W. I. Goldburg: Phys. Rev. Lett. 80, 277 (1998). 

[14] B. K. Martin, X. L. Wu, W. I. Goldburg and M. A. Rutgers: Phys. Rev. Lett. 80, 3964 (1998). 

[15] M. Rutgers: Phys. Rev. Lett. 81, 2244 (1998). 

[16] M. Rivera, P. Vorobieff and R. E. Ecke: Phys. Rev. Lett. 81, 1417 (1998). 

[17] P. Vorobieff, M. Rivera and R. E. Ecke: Phys. Fluids 11, 2167 (1999). 

[18] M. Rivera and X. L. Wu: Phys. Fluids 14, 3098 (2002). 

[19] M. Rivera, W. B. Daniel, S. Y. Chen and R. E. Ecke: Phys. Rev. Lett. 90, 104502 (2003). 

[20] G. Eyink: Nonlinearity 14, 787 (2001). 

[21] M.C. Lopes Filho, A.L. Mazzucato and H.J. Nussenzveig Lopes: Arch. Rational Mech. Anal. 
179, 353 (2006). 

[22] J. C. McWilliams: J. Fluid Mech. 146, 21 (1984). 

[23] R. Benzi, G. Paladin, S. Patternello, P. Santangelo and A. Vulpiani: J. Phys. A: Math. Gen. 
19, 3771 (1986). 

[24] M. E. Brachet, M. Meneguzzi and P. L. Sulem: Phys. Rev. Lett. 57, 683 (1986). 



19 



[25] M. E. Brachet, M. Meneguzzi, M. Politano and P. L. Sulem: J. Fluid Mech. 194, 333 (1988). 

[26] S. Kida: J. Phys. Soc. Japan 54, 2840 (1985). 

[27] K. Ohkitani: Phys. Fluids A 3, 1598 (1991). 

[28] K. Schneider and M. Farge: Phys. Rev. Lett. 95, 244502 (2005). 

[29] C. Basdevant, B. Legras, R. Sadourny and M. Beland: J. Atmos. Sci. 38, 2305 (1981). 
[30] B. Legras, P. Santangelo and R. Benzi: Europhys. Lett. 5, 37 (1988). 

[31] Y. K. Tsang, E. Ott, T. M. Antonsen and P. N. Guzdar: Phys. Rev. E 71, 066313 (2005). 

[32] M. Farge and M. Holschneider: Europhys. Lett. 15, 737 (1991). 

[33] J. Pedlosky: Geophysical Fluid Dynamics, II Ed., Springer Verlag (1987). 

[34] G. Paladin and A. Vulpiani: Phys. Rep. 156, 147 (1987). 

[35] Y. Jun and X. L. Wu: Phys. Rev. E 72 035302 (2002). 

[36] W.B. Daniel and M.A. Rutgers: Phys. Rev. Lett. 89, 134502 (2002). 

[37] H. Mizutani and T. Nakano: J. Phys. Soc. Japan 58, 1595 (1989). 

[38] R. Benzi and R. Scardovelli: Europhys. Lett. 29, 371 (1995). 

[39] B. K. Shivamoggi: Ann. Phys. 270, 263 (1998). 

[40] J. Weiss: Physica D 48, 273 (1991). 

[41] G. Falkovich and V. Lebedev: Phys. Rev. E 49, 1800 (1994). 

[42] K. Nam, E. Ott, T. M. Antonsen and P. N. Guzdar: Phys. Rev. Lett. 84, 5134 (2000). 

[43] P. G. Saffman: Studies Appl. Math. 50, 377 (1971). 

[44] B.K. Shivamoggi: Phys. Fluids 29, 769 (1986). 

[45] B.K. Shivamoggi: Phys. Lett. A 258, 131 (1999). 

[46] D.K. Rollins and B.K. Shivamoggi: Phys. Lett. A In Press (2007). 

[47] H. Kellay, C. H. Bruneau and X. L. Wu: Phys. Rev. Lett. 74, 3975 (1995). 

[48] J. R. Herring and J. C. McWilliams: J. Fluid Mech. 153, 229 (1985). 

[49] V. Borue: Phys. Rev. Lett. 71, 3967 (1993). 

[50] M. E. Maltrud and G. K. Vallis: J. Fluid Mech. 228, 321 (1991). 

[51] K. G. Oetzel and G. K. Vallis: Phys. Fluids 9, 2991 (1997). 

[52] M. Gharib and P. Derango: Physica 37D, 406 (1989). 

[53] H. Kellay, X. L. Wu and W. Goldburg: Phys. Rev. Lett. 74, 3975 (1995). 



20 



[54] J. Paret, M. C. Jullien and P. Tabeling: Phys. Rev. Lett. 83, 3418 (1999). 
[55] B. Dubrulle: Phys. Rev. Lett. 74, 959 (1994). 

[56] R. Benzi and L. Biferale: in Theories of Turbulence, Ed. M. Oberlack and F. Busse, Springer 
(2002). 

[57] A. N. Kolmogorov: Dokl. Akad. Nauk SSSR 32, 16 (1941). 
[58] R. H. Kraichnan: J. Fluid Mech. 62, 305 (1974). 
[59] A. N. Kolmogorov: J. Fluid Mech. 13, 82 (1962). 
[60] A. M. Obukhov: J. Fluid Mech. 13, 77 (1962). 

[61] U. Frisch, P. L. Sulem and M. Nelkin: J. Fluid Mech. 87, 719 (1978). 
[62] B. Castaing, Y. Gagne and E. J. Hopfinger: Physica D 46, 177 (1990). 
[63] K. R. Sreenivasan and C. Meneveau: Phys. Rev. A 38, 6287 (1988). 
[64] M. Nelkin: Phys. Rev. A 42, 7226 (1990). 

[65] R. Benzi, L. Biferale, G. Paladin, A. Vulpiani and M. Vergassola: Phys. Rev. Lett. 67, 2299 
(1991). 

[66] U. Frisch and Z. S. She: Fluid Dyn. Res. 8, 139 (1991). 

[67] G. L. Eyink: Physica D 91, 97 (1996). 

[68] E. Lindborg: J. Fluid Mech. 388, 259 (1999). 

[69] B. K. Shivamoggi: Ann. Phys. 312, 268 (2004). 

[70] G. L. Eyink: Phys. Rev. Lett. 74, 3800 (1995). 

[71] P. L. Sulem and U. Frisch: J. Fluid Mech. 72, 417 (1975). 

[72] A. Pouquet: J. Fluid Mech. 88, 1 (1978). 

[73] A. D. Gilbert: J. Fluid Mech. 193, 475 (1988). 

[74] B. K. Shivamoggi: J. Phys. A. Math. Gen. 23, 1689 (1990). 

[75] A.M. Polyakov: Nucl. Phys. B 396, 367 (1993). 

[76] P. Dmitruk and D.C. Montgomery: Phys. Fluids 17, 035114 (2005). 
[77] L. Kadanoff: Physica A 163, 1 (1990). 

[78] B. K. Shivamoggi: Chaos, Solitons and Fractals, 32, 628 (2007). 



21 




Figure la. Intermittency correction Aa*(p) vs. p for the intermittency parameter 
a = 10 (— 3D, - - - 2D). 
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Figure lb. Intermittency correction Aa»(p) vs. p for the intermittency parameter 
a = 100 (— 3D, - - - 2D). 
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